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1 Introduction 



AdS/CFT correspondence [1] relates conformal A/'=4 supersymmetric Yang-Mills theory 
(CFT) with string theory in AdS^ x space-time. Large number of colors N ^ oo and 
't Hooft coupling A = Qym^ ~^ ^ further lead to the classical supergravity regime (weak 
coupling) for the latter, putting the CFT into a strong coupling regime. In the decade 
since its invention, this correspondence became an indispensable theoretical tool, providing 
multiple interesting results about a strongly coupled regime of A/'=4 supersymmetric gauge 
theory. 

Among the earliest were calculation of the energy of a static electric dipole [2], based 
on a shape of "pending string" held at the AdS boundary at the positions of two static 
fundamental quarks, separated by distance L. For further reference we will need the EOM 
of the string'^ -'^ 



where Zm is the maximal string extension into z direction. We will also use notation L/2 = 
The resulting potential 



has the famous factor \/A (instead of A in the weak coupling Coulomb law. Soon this 
calculation was extended to include magnetic objects (monopoles and dyons) by Minahan 
[3], which can be viewed as an endpoints of the appropriate Di branes on the boundary. 
We will continue to discuss puzzles related with the dipoles in the next subsection. 

Naively one may interpret this answer by thinking of a strongly coupled vacuum as a 
dielectric medium, with a dielectric constant given by the ratio of the strong coupling result 
to the zero order Coulomb*^: potential 



Vcoalo,nU _ V^r(l/4)4 



.636\/A (3) 



^Maldacena StT'^ 

Although in a very qualitative sense this idea is not wrong, it is certainly not literally true. 
A proof of that are the calculations to be reported below, which shows that the stress tensor 



*'^Which is not the second order equation coming from the Lagrangian but the first (energy) integral of 

it. The eqn(l) represents half of the string, the other half is obtained by reflection x —x. 

*^Wc remind the reader that we include in it exchange dtio to scalars. It is equal to that from gauge field 

exchange for quark-antiquark pair, while for two quarks they have the opposite signs and cancel out. 
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distribution in space is very different from that in weak coupling. Of course, this is to be 
expected, as the strongly coupled vacuum recieves nonperturbative modification from the 
fields, leading to a nonlinear response. 

Now, a decade later, there is a spike of activity of using AdS/CFT to understand 
properties of the deconfined phase of QCD, known as Quark-Gluon Plasma (QGP) [4]. A 
number of phenomcnological considerations lead to a conjecture [5] that QGP to be in a 
'strongly coupled' regime (sQGP) at temperatures not too high above the deconfincmcnt 
temperature T = (1 — 2)Tc. It is in this domain where RHIC experiments at Brookhaven 
found a "perfect liquid" properties of sQGP. Among AdS/CFT-based works devoted to it 
are calculations of the energy loss [6] and stress tensor imprint [7] of the moving objects 
in thermal CFT plasma. Those are quite spectacular, providing in particular a compete 
picture justifying another hydrodynamical phenomenon, a "conical flow" in Mach direction 
around the jet. 

Lattice studies of sQGP have also indicated features indicative of a strong coupling 
regime. Those most relevant for this work obviously are studies of static charge pairs 
(electric or magnetic). Large deviation from a perturbative picture of a screened Coulom- 
bic potential are observed at T above the deconfinement transition Tc. More specifically, 
many features at T > l.STc suggest a "quasi-conformal" regime, in which all dimensional 
quantities (e.g. normalized energy density e/T^) show weak T-dependence. 

(Even larger deviations from perturbative approach - the Debye-screened charges 
- are seen for static dipoles at Tc < T < l.bTc. Here the entropy and potential energy 
associated with the string has very large part, linearly growing with distance in some range, 
see e.g. [10]. The question whether flux tubes - remnants of confining strings - can continue 
to exist in a plasma phase was recently studied in [11].) 

More generally, the dynamics of electric and magnetic gauge fields in a strongly cou- 
pled plasma remains very poorly understood. In particular, it has been suggested that 
sQGP contains large component of magnetically charged quasiparticles - monopoles and 
dyons, see [12, 13]. Studies of the energy distribution in plasma induced by static dipoles 
have been extensively done at zero temperature, demonstrating existence of quantum con- 
fining string: unfortunately similar calculations at T > Tc are not yet available. One may 
wander whether the deconfined QCD-like theories in those regimes are or are not similar to 
the vacuum of A/'=4 Gauge Theory at strong coupling. 

These ideas motivated our present calculation, in which we calculate stress tensor 
"imprint" of static dipoles in AdS/CFT. A simple diagrammatic picture of what is calcu- 
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lated is provide by Fig.2. Apart of the solutions themselves, to be given below in different 
regimes, there are few particular issues which would like to investigate: 

(i) How the dipole is seen at large distances r oo? What is the power of distance and its 

angular distribution? Can it be related to expected behavior of electric and scalar fields? 

(ii) What is the field near one of the charges? Can a non-singular part corresponding to 
the fields of a second charge and polarization cloud be identified? 

(iii) Is there a visible remnant of the Matsubara string, or a picture rather is of two polar- 
ization clouds? In particular, what is the r.m.s. transverse size ^J <y']_> at \y\ = (the 
middle point)? 

(iv) We will also consider an electric-magnetic pair: our main interest in that is to see if 
there is some nontrivial features related with electric-magnetic field interaction. 

1.1 Strongly coupled versus weakly coupled dipoles 

One issue discussed in literature (after AdS/CFT potentials been calculated) was whether 
some kind of diagram resummation can get the reduction*^ of the coefficient, from ~ A to 
~ \f\. Semenoff and Zarembo[14] have found that one can do so using ladder diagrams*^. 
Shuryak and Zahed [15] have noticed that such ladder diagrams in a strongly coupled regime 
imply a very short correlation time between colors of both charges 



which will be crucial for understanding of the large-distance field below. 

The "imprint" of the pending string on the boundary was first addressed by Callan 
and Guijosa [16], who had calculated an "image" due to scalar (dilaton) field propagating in 
the bulk. The boundary operator associated with a dilaton is trF"^. Our work is very close 
to theirs, except that we calculate much more cumbersome graviton propagation instead of 
a scalar one, to get the boundary stress tensor. 

Their main results was a distribution of scalar density at large distances from a dipole 
r S> L (i) has the form 



and (ii) is spherically symmetric. Both are very different from what one finds for the 
shape of the electric field of a weakly coupled electric dipole, which has (i) power 6 and 

*^We remind the reader that we discuss A ^ 1 regime. 

*''Which is exact for a round Wilson loop, approximate for rectangular ones. 



5t ~ L/A^^ 



(4) 



trF^{r) ~ L^/r^ 



(5) 
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(ii) has a characteristic dipole energy distribution {3cos'^{9) + 1) (9 is polar angle from a 
dipole direction). Our calculation to be reported also will show power 7 but will have more 
complicated angular distribution. 

The reason why the power is 7 rather than 6 was explained by Klebanov,Maldacena 
and Thorn [17]. Imagine Euclidean time and perturbative diagram, in which perturbative 
field of each charge can be written as a time integral over a propagator, from a world line 
of a charge to an observation point: it produces power 6. The nontrivial point is that in 
strongly coupled regime color time correlation [15] mentioned above require both charges to 
emit quanta at the same time; this changes a double time integral into a single one, adding 
one more power of the distance. 



2 Solving the linearized Einstein equations in AdS^ 

As is clear from Introduction, the source of gravity in our problem are strings extended into 
the AdS space. Naturally those are considered to be weak sources, so we will linearize the 
Einstein equations (with A = 6) 

Rixiy - ^RQhv + ^Qixu = -K^T^v (6) 

(6) and solve for small deviations from the unperturbed AdS^ metric. 
We choose to start with another form of (6): 

Rt^u +\k -——g^A = -n^T^, (7) 

with d = 5 Linearizion of the above gives: 

SR/j,^ - 4%^ = SS/j,^ (8) 

where 6S^^ = -k^ {6T^^ - ^Qixu) 

We denote weak gravity perturbation as dg^i, = h^i, and use an axial gauge in which 
the following components vanish hz^i = {fi = z,t ). We use the usual Poincare 

coordinates for the AdS metric: 



, 2 -dt'^ + dx'^ + dz"^ 

ds^ = (9) 



and set the AdS radius Lads to 1#^. 
*^Factors of Lads can be easily reinstated by dimensional analysis. 
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Expressing the modifications of curvature SR^i, in terms of h^i, , (see Appendix. A for 
a brief derivation) we get the following equations 



^h^z,z - ^hz = SSzz (10) 

2 {h,m ~ hra) = SSzm (H) 
1 Z 1 1 

— Cl/tj^jj + —hfnn^z 2 (^"^i" ~^ ^n,m) "I" ^ i.^,m,n l^mn^j^) ~ ^^mn (12) 

where we have defined h = g^"h\a, hm = g^'^h\m,a, ^ = ^"^ {~dt + d'i + d1) , and from now 
on latin indices stand for 4 boundary coordinates (m, n = t, x^, x^, x'^). 

We could in principle solve for h from (10), the result of which can help to solve for 
hm from (11). Finally solve for h„in with h,hm plugged in (12). However, we choose to do 
it in a slightly different way: As (11) is first order in z, it is only a constraint equation. 
With the boundary condition: hmn = 0(thus h = 0,hm = 0) at z = 0, we obtain 

hm = h^m - 2 / 6Szmdz (13) 
Jo 

(10) is second order in z, but it gives also a constraint when combined with (12): Denoting 
(m,n) as the mn component of (12), — {t,t) + Sj (a;*,a;*) gives: 
17 f 

-jh^z,z - -j^h^z = -SStt + ^^SS^^^^ -2 {-dSzt,t + SiiJ^^a-.^a..) dz (14) 
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Combining (10) and (14), we obtain the solution for h 

J dz- z (^SSzz + SStt - ^iSS^i^i + 2 J dz {-6Szt,t + ^iSS^x^xi)^ 



(15) 



With h obtained from (15) and hm eliminated, (12) becomes a closed eqn for remain- 
ing components: 

1 z 

— Dhmn "^hmn "I" '^h.mn,z ~ ^mn (1^) 

where a "generalized source" is Smn = SSmn - Jq {SSzm,n + SSzn,m) dz + ,m,n + 2^"mh,z 
The source terms created by the string are obtained from the Nambu-Goto action of 
the string in a standard way 

Sj^G = ^ / d^a I d^xJ-detgS^^^ (x-X (a)) 

27ra' J J 

2SSng 



ST"" 



-GSGf^u 
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here we use G^j^i, and gap to denote AdS metric and induced metric respectively. 
The string world sheet can be described by 



= X {t, z) , x} 







(18) 



The resulting source is as follows (we use the order of coordinate indices in the 
following 5-d matrices ^ and all absent entries are zeros) 

1 



55, 



fJ,U 



2TXOt 



j8 [x^ - x) 6 (x^) 6 (x3) 



Vl + a;2 - x'i 



/ 2x1+xl+l 
' 3 

XtXz 

-xt 



V 



XfXz 
3 

Xz 



-Xt 

-Xz 



-§ {xt 



XZ 



1) 



1 {^t 



1)/ 



(19) 



With (16), (15) and (19), we can solve for hmn, provided any explicit profile of the 
string. We will do this for three different string profiles separately in the following sections, 
and extract the corresponding stress tensors. 



3 The stress tensor of a static quark 



As a warm up, we will start with the case of a straight string, which corresponds to a single 
quark in A/'=4 SYM. The string profile is simply x {t, z) = 0. Substitute in (19), we obtain: 



V 



(20) 



3/ 



Static source leads to the metric perturbation hmn which is time-independent. Per- 
forming a Fourier transform = J hmn^^'^^d^x we convert the PDE (16) to an ODE: 



mn.z "mn 



(21) 



An upper index k will be used below to indicate a Fourier transformed quantity. 
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S^jj is just S^i, without delta functions, hf' and s^,^ have simple forms displayed as 



follows: 



9 27ra' 



' 3 



V 



Z + 



J 



/o 



V 



J 



(22) 



(23) 



The equation is Bessel type and can be dealt with using a Green function built out of such 
functions. Instead we consider a more general equation with arbitrary power of z in the 
source 



-^2 ( 

2 \'' 



""mn J 



Oh'' -I- 

mn ' mn,z 



CnZ 



(24) 



which is directly solvable in terms of Meijer-G funcion and hyper geometric function: 

')^2^2 1 \ 2^-1 ■ 



+K2{kz) (Ci -1 F2 



f+l,f-l,0 



K2,3 



k^Z^ 



An + 8 



(25) 



The constants Ci and C2 are to be fixed by boundary conditions. One of the condition 
is the metric perturbation vanishes at AdS boundary, i.e. = at z = 0, which fixes 

Ci = 0. The other boundary condition proposed in [9] for thermal AdS is incoming metric 
perturbation at the horizon. However in our case, we need a different boundary condition 
due to the absence of horizon in AdS. Since /i^ grows as z^ in the present case, while 
show possible exponential growth at large z. It is natural to propose no exponential growth 
at z = 00 as the boundary condition. 

At large z, only the first term containing l2{kz) is dominant, the boundary condition 
becomes: 

m-i _ , / 1,2^2 1 



r< I r^2,i / k z 







(26) 



The asymptotic of Meijer-G function {z 

'k^zh^ 



-^2,1 
^^1,3 



4 |2+i,^_i,o 



00) gives: 



(27) 
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which finally fixes C2 = — ^^^^^T (| + l)r(^-l). Applying it to our source = 
ciz + C3Z'\ where ci and C3 are matrix- valued (the indices are suppressed here), we have 

TT „ 37r „ 

The stress tensor of the corresponding boundary CFT is proportional to the coefficient 
of z"^ term, which we denote as Qmn throughout this paper, in small z expansion of hmn- 
The precise relation can be obtained from (36) of [7], which in our case is simply (with 
Lads = 1): 

2 



Note Glj ^ 



-+i --1 > 

2^ '2 ' 



1,2 2 \ 

and 1F2 yn_^2 3^^~J -^"^^^ contains only odd power of z 
for odd n, thus does not contribute to Qmn- We have 

Qmn = lk^C2 (29) 

Reinstate the factor L\j^g, together with the relation ^^^^ = \/A, we have the final 
stress tensor: 



TT 



\ 



TT 

k 



V 



kni^n 



TT 



(30) 



It is easy to verify the stress tensor above is traceless T^nV"^"' = 0) which is a conse- 
quence of conformal invariance. It also satisfies the conservation of energy and momentum 
k"^Tmn = 0. In doing inverse Fourier transform, we find the k-integrals are not well-defined. 
One trick is to introduce a regulator e""*' (a > 0) to the integral, and take the limit a ^ 
in the final answer. We end up with the following result: 



-Vx 1 

TT StT 



2 
'3^ 



/l 



V 



+ 



V 



/o 



V 



yrnVn 



4 

3r6 



(31) 



*®We remind the reader that unperturbed metric has l/z^ and thus the relative smallness is 0{z^) fitting 
the dimension of the stress tensor. 
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where r is the distance from the quark. The ^ power is obvious by dimension. Let us recall 
the result obtained in [16]*'' 



1 



327r2 r4 

While in our case, the Too component gives 

1 



(32) 



1 



(33) 



127r2 

In both (32) and (33), the dots represent contributions from scalars and gluinoes. 
If we assume the magnetic field is not present, the difference in the two operators implies 
significant contribution are received from the scalars and gluinoes. 



4 The stress tensor image of static electric dipole 

Now we turn to the Maldacena's pending string, the ends of which attached to a quark and 
antiquark, corresponding to a static electric dipolc. The string profile x{z) is double- valued. 
We use ±a;(2;) {x{z) > 0) to denote two halves of the string. The EOM can be integrated 
to give x{z) in terms of elliptic integrals. We will not refer to explicit form until the end of 
the calculation. 

The source term and its Fourier transformed version are a bit complicated: 

1 



dS, 



— K^Z 



3 



3 

2^2 



j5 [x^ - x{z)) 5 (a;2) 6 {x^) 



\ 



Xz 

-xl+2 



2x1+2 



2xi 



+ {x 



-X) 



(34) 



*'^There is a typo in eqn (23) of the paper. We quote the corrected expression 
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dS 



27ra' ^ 



2x1-1 



2x2+2 



2xi+2 



cos(A:ix) 













\ 












Xz 






+ 





Xz 







i sin(A;ix) 

















I 






Oy 





(35) 



It is understood that the source term vanishes for z > Zm- (15) and (16) gives: 



h{z < Zn 
h{z > z„ 



1 -K- 



3 27ra' 

1 -K^ 



3 27ra' 



7{nzm) + ,:iz'-zi)Gizr, 



(36) 



/l \ 



27ra' 



-1 



/l \ 



+ 3£^2(^) 



V y 



+H{z) 



/o \ 

A;2 

V ^3 y 



2 



/o 



27ra' 



kmkn 



G{z„ 



\ J 



2ki k2 fca 
A;2 

V ^3 

/-I \ 



V 1/ 



-1 



V 2/ 



/O \ 



+ 



F'jz) 
6z 



V / 



/-I \ 



V V 



(37) 



F{z„,) + -iz^-zl)G{zm) 



(38) 
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with 



E\{z) = x\ cos(fcia;) 



2z 

E2{Z) = r—^ COs(fcix) 
V + 



-4cos(A;ia:) n , , f — 4sin(A;ia;) , 



2 



V 



Jo ^/TT^, Jo \ Jo ^/T+ 

G{zm) = I -^^^t^zkiXzdz 

Jo yJl + X% 

2z sin(A;ix 

-X 

With the exphcit expression of we can build the general solution to (16): 

hi, = h{kz)C2 + K2{kz)C^ + 2 (h{kz) r id^&^dz - K2{kz) r 'IrMMh^ldz) 

At large z, no exponential growth condition requires C2 + 2 imnMI^^iMldz = 0. 
The convergence of the integral is ensured by K2{kz) in the integrand . At small z, ~ 
0{z) while hikz) ~ 0{z^) the integral containing hikz) is finite as z approach 0, therefore 
the boundary condition gives: Ci - 2 fLiMMMlcf^; = 0. 

In order to extract the stress tensor, we need to collect terms. It is helpful to write 
down the series expansion of the two integrals 



J^z, 



'iiiMsMrf, = + + (41) 

= ^ + ;,, + ... (42) 

z Z Z 

t2 , 



The coefficient of z^ is given by: Qmn = ^^^^2 + ^bo. Note Ci does not appear in 
the expression We may also write it as 

^ J Zm ^ ^ \Jzm ^ ^ J 

We could proceed in momentum space. However it turns out to be much easier and 
illustrating to do inverse Fourier transform and continue in configuration space from now 
on. 

A nice property of Fourier transform is T~^{F{k)G{k)) = J f{x)g{y — x)dx. Identi- 
fying the source dependent as F{k), the inverse Fourier transform of which gives f{x). 
Correspondingly, each G{k) is transformed to g{y — x). The latter can be interpreted as a 
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propagator from a point on the source x to a point on the boundary y. With this in mind, 
we define the following propagator: 



Ps{y- ^) = 7A3 / k^K2{kz)e-'^^d^k = ^ ^ (44) 

Let us take a moment to worry about the term involving h-\. By analyzing small z 
behavior of Smn and l2{kz), we find b-i = Inverse Fourier transform of ^6-1 is not 
well-defined. Again we introduce the same regulator e~°'^ as in the previous section. We 
find a vanishing result after taking the limit a — >■ 

Finally, we can write the stress tensor in a very short form: 

+ \ f dz [ ^mn{z,x)Ps{y-x) ^,^ ^^^^ 

Before proceeding with the calculation, we would like to make few general comments: 
(i) The trace of the stress tensor is given by the coefficient of term of h. From (36), we 
find that h ^ F{z) at small z does not contain 2;^ term, therefore we expect the final stress 
tensor to be traceless, which is also required by conformal invariance. (ii) The divergence 
of the stress tensor d\Txjn turns out to be the the coefficient of term of hm- From 
(13) and (36) we conclude dxTxm = {5{x^ - Xm) - + Xm)) 5{x^)5{x^)5mi. The 

divergence is non- vanishing only for m = at the end points of the string where the quark 
and antiquark are placed. It corresponds to a pair of antiparallel forces which hold quark 
and antiquark, preventing them from falling onto each other. This will be another general 
condition to be satisfied by the stress tensor. 

4.1 Far field 

With (45) at hand, wc first calculate the stress tensor in region far from the dipole. The 
inverse Fourier transforms of s„j„ are linear combinations of those of Ei,E2, F, G, H. Such 
terms as kmH can be replaced by idx^H = —id^mH = idy^. In the first identity, we 
use partial integration so that the derivative only acts on the propagators (we indicate this 
with a left arrow on top of the derivative) . The second identity is due to Pg = Ps{y — x). 
Similarly, kmknF -"d^Jd^^F -'dy^'dy^F 



**this may seems problematic. Actually the same regularization can also be applied to K2{kz) if we first 
expand it in series of k. The non- vanishing terms match those obtained from series expansion of propagator 
P, in r 
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We list the back-transformed result of Ei, E2,F, G, H here: 

v5 



El 

E2 = 
F = - 
G 



z 



-.S{x^)^{x^)^{x^ - x{^z)) + {x^ ^ -x^) 



-5(x2)<5(x3)(5(x^ - x{z)) + {x^ ^ -x^) 



— b{x^)b{3? 

4rS(x'^)5(x^ 



-(z^-zf){zi-3z 



4zl 



^0{z-zi] 



(x 



-x 



2z\ 



^'ieiz-zi) + 



4 _ 4 



2zi 



^6{z-zi) +{x^ 



H = -^Six'')Six^ f''^'^ '' 9{z - zi) - ix' 



(46) 

(47) 

(48) 
-x^) (49) 
(50) 



with zi = zi{x^) (0 < < Xm)-, the inverse function of x{z). Contribution from 
negative x^ is included in the second term for each function. Note Ei,E2,F,G are symmetric 
in x^, while H is antisymmetric. 

In order to obtain the far field stress tensor, we need to perform a large \y\ expansion 
of the stress tensor. Note the y-dependence enters the stress tensor via the propagator, we 
can do a large \y\ expansion on the propagator in the second term since z < Zm <^ \y\- 
While for the first integral, z extending to infinity, we need to do the integral first before a 
valid expansion is possible. Fortunately this time the source has very simple z-dependence: 
Sjnn{z) = # + #z^. The rest of the calculation is straight forward. After collecting all 
terms, we find the first nontrivial result appears at the order The power again agree 
with the result of trF^ obtained in [16]. We list the stress tensor as follows (up to the order 



^ 1 -\/A 15 / aG{7yi - y'^) /aEi aE2 Qf\ 1 
°° ~ 4 TT 47r V 12|y|9 V 3 3 6 / 



'■Om 







1-VA 15 
4 TT 47r 



-7aH + -ac 
6 



'^2yi 2/2 2/3 

2/2 

\y3 



/l o\ 



(51) 

(52) 



^2 1 

-^ + ^{aEl + aE2)dmn-^ 



1 (Tap 7aG\ ymVu '^'i-ac ymynVi 



21 



\y\' V 6 12 J 4 \y\^^ 



(53) 
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with 



aG = 2 G{zm) {x^) dx^ = -0.7189^; 
Jo 

{"■Cm 

aF = 2 F{zm)dx^ = -0.95852;^ 
Jo 

aH = 2 H{zm)x^dx^ = 0.1797z^ 
Jo 

fO 

asi = / Ei{z)dz = -0.7189z^ 

J Zm 

fO 

aE2 = / E2{z)dz = -0.47934 

" Zm. 



(54) 

We can verifiy explicitly that the stress tensor is traceless and divergence-free at this 

order. 

Now we proceed to analysis of the results, describing which features are general and 
should be expected and which of them are qualitatively new. 

A vanishing energy flux (Poynting vector) Tq^ = is related with zero magnetic field 
expected for static electric configuration. Indeed, a time reversal would change the sign of 
the magnetic field and the Poynting vector, but leaves the problem invariant. 

Having said that, we by no means imply that the only field in question is the elec- 
tric field. Indeed, vacuum polarization should include all other fields of the theory, and 
perturbatively we know that all color fields of the theory - gluinoes and scalars - should 
contribute, to charge polarization density as well as to the energy we calculate. However, a 
very simplistic view of the scalars*^ based on (d^^)^ Lagrangian would produce the same 
distributions as a vector field, since that can be viewed as just generated by another scalar 
field ^0- 

The obvious point of comparison is stress tensor distribution for a perturbative dipole. 
Its electric field 

F ( yrn-{L/2)em yrn + {L/2)e^ \ 

^^7:>Wr,-{L/2)em? \ym + {L/2)e^\^) ^''''> 

leads to stress tensor which is at large distances ^ I? j'lf' . The result we obtain is ~ jr^ : 
the difference is due to the a phenomenon of "short-time-color-locking" [15, 17] we already 
discussed in the Introduction. Perhaps another way to explain it is to say that a scalar 
density, induced by a dipole, is large in all the volume ~ L^. 



*®Ignoring quartic terms with commutators of various flavor components. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 



theta 

Figure 1: (Color online) The far field energy distribution in polar angle 6{cos{9) = yi/\y\), 
normalized at zero angle. Solid (black) line is our result, compared to the perturbative 
result (Scos^ + l)/4 given by the dashed (blue) line. 

Let us now comment on the angular distribution. Perturbative dipole field at large 
distances contains the first power of the dipole vector: thus its angular momentum is 1. 
Energy density constructed out of this field, obviously has only angular momenta 2 and 
0, or powers of cos"(^) with n = 0,2. Stress tensor also contains such components, but 
also terms of the type Tmn ~ ymUniLyY ■ Looking at our result we find that indeed no 
other angular structures appeared. This is to be expected, as electric field is still the only 
vector field of the theory. The angular distribution of the far field energy is compared to 
the perturbative result in Fig.l: although there is tendency to a more spherical distribution 
(like obtained for scalar density [16]), the peaks in the dipole directions are still there. 

One more simple case to discuss is the stress tensor on a line connecting the charges: 
by symmetry transverse component of the field E±_ = and only remains. The 
Maxwellian tensor then should satisfy T22 = = — Tn = Tqo: and the result we ob- 
tain does not satisfy it. We thus see once again, that gluino and scalar parts of the stress 
tensor must contribute to the far field asymptotic in question. 
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4.2 A field near one charge 

Next we would like to study the stress tensor near one of the charge. For this purpose, we 
make a shift of variables yi — ^ yi + Xm, y2 — ^ 2/2, 2/3 — ^ 2/3, and consider small \y\ behavior 
of the stress tensor. 

It is clear from the single charge result(31) the stress tensor will blow up as \y\ — 0, 
so to the leading order in \y\, we may focus on its divergent part only. 
Let us recall the basic expression for the stress tensor: 

Qmn = -7 / dz ^ ^ -d-'x 

+ \[dzj 'mn{z,x)Ps{y-x) ^,^ ^^g^ 

As \y\ — 0, the first term is finite(2; > Zm), which we ignore as discussed above. While 
the propagator in the second term Pg = — — — r contains a singularity at z = 0, r = 0, 
which leads to a possible divergence in stress tensor(unless the source provide enough powers 
of z). We can also claim the divergence is from integration at small z. Since the integral 
involving Smn and Pg cannot be done analytically, a careful analysis is needed to obtain the 
leading terms in Laurent expansion of the stress tensor. 

We first use the common factor 5{x'^)5{x^) in the source to simplify the propagator: 

= (yi + xm- x^f + (2/2 - x'^f + (2/3 - x^f = (yi - ^xf + yl + yl 
= To - 2yiAx + Ax^ 

with = yl + y2 + y"^-, Ax = x^ — Xm- Then the propagator can be expanded in Aa;-^: 

= ^ f ^ - 7 ^yi Ax^ + • • • 1 (57) 



47r 



[z"^ + rl) 2 (2;2 + rl) 

Note the leading term of the propagator does not depend on x^, x^, x^. A similar trick 
is used as in the case of far field: km = —idxm = ''■^ym- The second identity is due to 
Pg = Ps{y — x). If only the leading order result of the stress tensor is needed, we perform 
the x-integral with the source, keeping the smallest power in z(As we argued before smaller 
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power of z corresponds to larger term in expansion of the stress tensor): 

^5 



/ 



Eid-'x j- 



X ^ z 



J Fd^x ~ -4 j {z^ - zl)B{z - zi)dx^ = -\ 

/F' z 
—d^x 2^-2 6'(z - zi)dx^ = -2z 

j -iHd^x ~ zid{z - zi)dx^ 



Convolute the above results with the leading order propagator, we find they give the 
following divergence: 

El -> ln{ro) 

^ F' 1 

E'2, > ^ 

z 

ro 

Therefore the leading order result is given by £^2, "7- and F. The last also give \ when 
combined with the double derivatives in the coefficient. Collecting all the contributions, we 
find the leading near field contribution, which is of course precisely the stress tensor of a 
single charge (31) 

/l \ /O \ 



rr^LO ^ 1 



34 



\ 



+ 



yrnVn 



V V 



3r6 



(58) 



The aim now is to extend the analysis to the next order correction to (58). Note the 
correction from the source will give at least O(z^) correction, while that from the propagator 



is of 0{Ax^] 



~ Zi ~ 



z^, with an additional + Tq in the denominator. As a result, we 



can keep the leading order source but care about the correction from the propagator when 
necessary. Finally we find the next order correction to the stress tensor is from the LO 



source E2,F,^ convoluted with the NLO correction to the propagator —7 
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as well as the leading result from iH. We display the correction to the near field as follows: 

/5 \ /O \ 



-,NLO 



VA J 1_ 



8/ 



2yi y2 ys 

Z/3 



3rl 



J 



/o 



V 



yrnVn 



yi 



(59) 



We can also verify the stress tensor at this order is traceless and divergence-free. 

Let us now analyze the results and compare it with expectations. In general one can 
expect that close to the charge there is a singular electric field E^'^'^^ ~ l/rg plus a finite 
field induced by all other charges. 



EiEj 



+ 



(60) 



The scalar field in weak coupling add the same distributions. 

If the vacuum would be a simple dielectric, both the singular and regular field would 
be just free fields times the dielectric constant (3), and the relative correction be the same. 
Let us see whether this idea works or not. In weak coupling^^'' the correction to Tqq is 
1 -2(?/ir)/L2 while our strong coupling result gives 

,2(yir) 



Tqq ^ _ {yir) 

rpLO y2 

-'OO m 



1 - 0.34- 



L2 



(61) 



The sign and the structure of the local field is the same, while the magnitude is addition- 
ally reduced by about a factor 1/3. What we learn from this comparison, once again, is 
that although a strongly coupled vacuum of the theory works as a polarizable dielectric 
qualitatively, this is not true literally. 



4.3 Is there a visible trace of the string? 

Another interesting question is the transverse distribution of energy. In particular we calcu- 



late the r.m.s.: y^< y| > = ^ '^j^Too'dyT ^ =ya=oj " ' which characterizes the transverse energy 
distribution on the middle plane between the quark-antiquark pair. 



#10 



There axe both gauge and scalar fields, but distributions they produced in zeroth order are the same. 
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1 f soo{z,x)Psiy- x) 



Too ~ Qoo = 4 / dz I ""^ ' ^^^^^ '-d'^x (62) 



, , Ei{z) E2{z) F'{z) 

.oo(.<..)~^ + ^-^ 

G{Zm) 



S00{z > Zm) 



6 



Note y-dependence enters only through the propagator Pg, we can do the y2 integral 
with the propagator first, then convolute the result with the source soo- The rest of the 



calculation is straight forward. We will skip the details and only give the result: \/< > ~ 
OAlzm, while half the size of the dipole is « O.GOzm- The r.m.s. is about ^ of the dipole 
size, smaller than the perturbative result \J <y\> = ^ ■ 

In order to make the trace of string clear, we would like to rewrite (45) in a more 
physical form. This is done by defining: s!^^ = Smn — Smn, then we have 



CXD 



, S^„(z,x)Ps{v — x) ,q , , 

dz / "^"^ ' ^ -d^x (63) 



^JO 

The first piece is sourced by the original string Smm while the second piece corre- 
sponds to contribution from sj^„. Since the latter is obtained from S^i, via (15) and (13). 
The transform from S^^ to can be interpreted as a bulk-to-bulk propagator, which is 
then attached to the bulk-to-boundary propagator Pg to contribute to the stress tensor. We 
schematically illustrate the two contributions in Fig.2 

We use the component Too as an example to study the relative contribution from the 
two pieces: 



6 Jo 



^4 _^4 



7 

[z^ -I- a;(z)^)2 



=_! r dz r"^ dx^ ^^~^^^ 9{z-zi 



6 Jo 'Jo "" zl (z2 + (xi)2)i 



-1 zi-3zf 1 



6 Jo V4^izlb{zl + x\ziff2 



(64) 



(65) 
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Figure 2: (color online) Schematic demonstration of the pending string and the propagators 
of stress tensor. The source is at the point A integrated over string, it either (a) goes 
directly to the observation point y via bulk-to-boundary propagates(dashed line) , or (b) 
first transforms to s!^^ in some other point B via bulk-to-bulk propagator (dash-dotted line), 
then goes to the observation point 

We plot the integrands of (64) in Fig.3. All three curves have a peak at z — Zj^i, 
which is due to geometry of the string. However the peaks are square root singularities of 
geometric origin, which do not contribute significantly to the integral and the finial Too- 
Instead the latter receives significant contribution from integration of all values of z. 

5 A field of electric-magnetic dipole 

It is also interesting to consider the stress tensor of an quark and monopole, in which case 
both electric and magnetic fields are obviously present. The string profile of the electric- 
magnetic dipole is obtained by Minahan[3]. It consists of a (1, 0) and a (0, 1) string, attached 
to the quark and monopole at z = respectively, and a (1, 1) string extending from z = oo. 
The three string attach to each other aX z = zq, forming a Y-junction. With a suitable 
choice of coordinate, we can describe the (1,1) string by = 0, and describe the (1,0) 
string and (0,1) string profile by = x{zmi,z) > and x^ = —x{zm2,z) < 0, where 

2 

ZTni,Zm2 are parameters of the string profile x{z). x{zm,z) satisfies Xz = , ^ The 

parameters given by [3] are: 
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Figure 3: (color online) The integrands of the z integral along the string for Q\q{\A\\& 
dotted), Qoo (gi'een dashed) and their sum(red solid), with = 1 



Zm\ = z^ax (66) 

Zm2 = 2:oa2 (67) 



ai=(^2— ) a2 = (l + i^)* 



where t = ^, g is the string coupling. 

The action of a (p, q) string is given by: 

+ ^2^2 



S^G = - ^^^^"^^f ^ / dVy^ (68) 
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The SSf^i, following from the action is: 



' 3 



— K^Z 



Xz{Zml,z)'^+'i 

3 



2Xz{Zml,zf+2 

3 



V 

+ {x{Zml, Z) -x{Zm2,z)) t 



2Xz{Zml,z)^+2 



(I 



-K^Z 



5S^,{Z > zo) = -^S {x^) S (x^) S {x') 



\ 



VT+t^ (69) 



We are not going to elaborate the calculation in any detail, as the same procedures 
for the electric dipole's case apply. The far-field answer is 



-\/AJ_ 

TT StT 



/l 



V 



\ 



1/ 



/o 



V 



Umyn 



\ 



J 



(70) 



Difi'erent from the case of electric dipole, in which the leading total charge term drops out, 
the leading order now is . The result is proportional to ^/\(T+W) = ^^{gYf^j + {AttY /g- 
in good agreement with electric-magnetic duality*^^ of the problem. This shows the electric- 
magnetic dipole looks like a dyon to distant observer. 

Perturbatively one expect no correlation between electric and magnetic charges, and 
the answer proportional to a sum A?^(5yj^^/47r-|-47r/g'yj\^), without a square root. The reason 
a common square root appears can again be traced to color correlation time by Shuryak 
and Zahed: for example they have also shown that Coulomb, spin-spin and spin-orbit forces 
are also united into one common square root [20]. 

For the near field, we recall the calculation of the previous section, the LO stress 
tensor near the quark(monopole) is again the same as that of a single quark(monopole) . the 



km)' 



*^^We remind the reader that Dirac condition in this theory is simply that magnetic charge is the inverse 
of the electric one. 



23 



NLO stress tensor only depend on the profile of the string attached to the quark (monopole). 
Therefore, we can obtain the stress tensor by the substitution: Zm ^^mi (quark), Zm —>■ 
Zm2 (monopole) . We display the NLO near field result for the quark and monopole in 
(71),(72). 



/o 



-Vx 1 ^_ 

TT 127r 4i 



yrnVn 



6rQ 



yi 



/5 \ /O \ 

2yi 2/2 ys 
y2 

sy V 2/3 J 



34 



(71) 



/5 



-Vx 1 t 

^ 127r ^2 



yi 

6rQ 



+ 



/o 



V 



yi 



I 



\ /o 



+ 



8/ 



22/1 2/2 2/3 
2/2 
2/3 



3r3 



(72) 



The result at NLO suggests the impact of a monopole to the quark is the same as 
an antiquark at some distance away. The precise relation between the quark-monopole 
distance Lqm and quark-antiquark distance Lqq can be estimated. Lqq should be chosen 
such that Zm reproduce z^i for Lqm- (2.6) and (3.2) of [3] gives: 



00 



dy 



Lqm = zo ^cti J 



1 2/Vy^ - 1 

°" dy 

ai 2/2 



+ ^2 



zoai 



i: 



dy 



2/V2/^-l 



dy 



0-2 



yWy^ - 1 



(73) 
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where the approximation is due to the Umit g 0, t —>■ oo. The result shows that 
in the above Umit, the quark feels the monopole like an antiquark at twice the distance. 
Similarly, the monopole feels the quark at a distance L like an antimonopole = -^L 
away. 

Finally, let us address the ussue of the angular momentum and Poynting vector. Pcr- 
turbativc chargc-monopole pair has at a generic point electric and magnetic fields crossing 
at some angle, thus producing a nonzero Poynting vector Tq^ 7^ 0. In fact its direction is 
rotating around the line connecting charges, leading to nonzero angular momentum of the 
field. In fact the Dirac quantization condition is known to be directly related to quantization 
of this angular momentum. 

However, in our setting with Minahan's solution this effect is entirely absent and 
there is no anular momentum or Poynting vector, = 0. This can be traced directly 
to the expression (70) for the source which has no such component. In gravity setting the 
energy-momentum of the Minahan string construction does not care about direction of the 
magnetic flux, and the problem is again static and t-reflection symmetric. 

Perhaps the way to remedy the situation is to start with a different classical rotating 
string, with some nonzero angular momentum, which value is to be tuned to fit the Dirac 
condition. If we will be able to make progress along this line, we will report it elswhere*^^. 



6 Summary and outlook 

The main results of this work are general expressions for the stress tensor induced by 
objects in the AdS bulk (31), (51), (58), (59), (70), (71), (72). In general, we found that two 
components of gravity perturbation - the trace of the metric h and its tensor part /t^,/ 
- have different equations and Green functions. Although h itself on the boundary does 
not have 0{z'^) corrections or induced stress tensor (as follows from conformal symmetry 
of the boundary theory), two components are intermixed in curved background and thus 
h (incorporated into a "generalized source") leads to physical effects including the stress 
tensor. 

General formulae are then used for static electric and electric-magnetic dipolcs, as 
important examples. Confidence in the results come from checking all of them for traceless- 
ness and energy-momentum conservation. We worked out the far field asymptotic, as well 
as an expressions for the field near one of the charges. 

'^^^We thank Andrei Parnachev and Jinfeng Liao for helpful discussions of this issue. 



25 



The far distance asymptotic of the stress tensor is ~ Li^ /r^ ■, the same as in previous 
calculation [16] for dim-4 scalar density, the angular distribution is different. We found 
that although all angular structures are as expected from perturbative analysis for dipoles, 

the coefficients (and angular distribution of stress tensor) arc quite different from the weak 
coupling limit. The same is found for the near-field domain. It means although a naive idea 
of strongly coupled vacuum acting as a dielectric qualitatively is holding, quantitatively it 
definitely fails. 

We also found that on the boundary there seems to be no visible trace of a string. 
In fact even in between the two charges (e.g. at yi = 0) the dominant contribution still 
comes from "vertical" parts of the string rather than its "horizontal" part directly beneath 
the observation point. The distribution looks like two distorted polarization clouds about 
two charges, instead of a string-like object. 

This conclusion is relevant for interpretation of string-like entity which seems to 
appear via linear part in static dipole potentials on the lattice at T just above deconfinement 
for QCD-like theories. We think those are due to some flux tubes. Their formation is due to 
phenomena which needs more specific ingredients than just a strong coupling regime. Let 
us further conjecture that the distributions we calculated from AdS/CFT should instead be 
similar to those in QCD-like theories in a "quasi-conformal regime", at temperatures not 
too close to deconfinement, T > 1.5Tc. This is the region in which flux tube effects are 
gone, the potentials become a screened-Coulomb type and thermodynamical observables 
are about constant when divided by appropriate powers of T. This conjecture will be tested 
directly in forthcoming lattice calculations. 

As an outlook for this work we have in mind, we would like to work out stress tensor 
imprints of dynamical (rather than static) objects. In particular, those arc "debris" created 
in high energy heavy ion collisions, see [18] for a basic picture and to our previous paper 
[19] in which we formulated the picture and calcTilatcd trajectories of different types of 
objects falling into AdS bulk. We hope then elucidate the process of black hole formation, 
out of those "debris" and see whether the stress tensor imprints would be approaching 
hydrodynamical solutions, which were so successful for the description [5] of RHIC data. 
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A Linearization of Ricci tensor 



We may start with the following relations: 

Since the covariant derivative on the metric vanishes, g^^^ = 0, the metric commutes 
with the covariant derivative. 6R^i, can be further simplified. 

^^ixv — 2^ ( hcrfji-^v,\ ^cri/;/i;A "I" ^crA;/i;i^ "I" ^/u;/;(t;a) 

= -\g^%Ki.;v,x + K.;i.-x) + ^h,^;. + ^5^"V;^;A (74) 

with hfj,^ = 5g^„ h = g^'^hxa- 

We choose to work in Poincare coordinate, the only nonvanishing Christoffels of which 

are: 



T^z _ -pz _ ^ -pz _ -'■ _ pa;' _ -'■ ('7c\ 

^tt—^ZZ— ^ xixi — ^tz — ^ xiz — ^ \'^) 



We calculate the components SRzz,SRzm,SRmn separately. Through tedious algebra, we 
arrive at: 

SRzz = ^hz,z - ^hz (76) 

^Rzm = - {h^m — hm)^z (77) 

1 z 1 1 

^^mn ~ Q^^^mn "I" "^h^pim + ^^mrijZ ^ (^m,n "I" ^n,m) "I" ^ (^,m,n ^mn^,z) (78) 



with hm = g^''hxm,cT 
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